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There has been an increasing interest in the monitoring of nuclear fuel for power reactors by 
detecting the anti-neutrinos produced during operation. Small liquid scintillator detectors have 
already demonstrated sensitivity to operational power levels, but more sensitive monitoring requires 
improvements in the efficiency and uniformity of these detectors. In this work, we use a montecarlo 
simulation to investigate the detector performance of four different detector configurations. Based 
on the analysis of neutron detection efficiency and positron energy response, we find that the optimal 
detector design will depend on the goals and restrictions of the specific installation or application. 
We have attempted to present the relevant information so that future detector development can 
proceed in a profitable direction. 
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I. INTRODUCTION 

Recently, there has been an increasing inter- 
est in the monitoring and safeguarding of nuclear 
power reactors. Over the last few years, a group 
from Livermore and Sandia National Laboratories 
has been demonstrating the feasibility of a small 
simple detector to monitor the anti-neutrino pro- 
duction during the operation of a nuclear power 
reactor located in San Onofre, California lj. By 
using the anti-neutrinos which are produced by the 
uranium and plutonium fuel itself during the fis- 
sion process, the reactor fuel can be continuously 
monitored in a non-invasive way. The ability to 
monitor the nuclear fuel composition in real-time 
has advantages for both limiting the proliferation 
of nuclear material as well as increasing the oper- 
ational efficiency of power generation[20- 

The current experience at San Onofre has 
demonstrated that a small neutrino detector lo- 
cated within 25 m of the reactor core is easily sen- 
sitive to the power level at which the nuclear re- 
actor is being operated. However, more sensitive 
tests to determine the fuel burn-up and fuel com- 
position, while showing promising results, have 
shown the limitations of the detector design. 

This work is an attempt to evaluate several pos- 
sible directions for improved compact neutrino de- 
tector designs using the latest simulations avail- 
able to the reactor neutrino community. The de- 
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tectors are based on the conventional technology of 
liquid scintillator and photomultiplier tubes. The 
anti-neutrino event signature is the inverse beta- 
decay process: D e + p — > e + + n. This yields 
a coincident event signature, from the prompt 
positron annihilation and the delayed neutron cap- 
ture, which is relatively free from background con- 
tamination. The use of a liquid scintillator that 
is doped with gadolinium improves the signal to 
background further by reducing the neutron cap- 
ture time and increasing the energy which that 
capture releases. 

The simulations were based on the open-source 
libraries of the Generic Liquid- scintillator 
anti-neutrino detector Geant4 simulation 
(GLG4sim) 4, 5] which have been extensively used 
in the KamLAND experiment [(|. In addition, we 
have made use of several of the improvements that 
have been developed within the Double Chooz 
collaboration^}. Specifically, the development 
of high quality Gd doped liquid scintillator has 
been carefully studied and implemented within 
the simulation, as well as detailed improvements 
in material properties and PMT performance. 

For this study, it was considered that there are 
two major features of the detector response which 
will be important for reactor monitoring. The 
first concerns the efficiency for detecting events 
in which the neutron is captured on Gd within 
the fiducial volume of the detector. The overall 
acceptance for anti-neutrino events is a product 
of positron and neutron identifications. However, 
the positron detection efficiency is largely a func- 
tion of the environmental conditions and can more 
directly be improved by shielding the detector to 



allow a reduction in the minimum energy thresh- 
old for acceptance. The neutron efficiency, on the 
other hand, can be improved through design fea- 
tures of the detector — ensuring that the photons 
released through the n-Gd capture are all detected 
within the scintillator to allow good separation be- 
tween neutron-like and positron-like energy depo- 
sitions. 

The second detector performance characteristic 
studied was the uniformity of energy response for 
positron events. If sufficient statistical power is 
available, the positron energy spectrum — which is 
directly related to the parent neutrino spectrum — 
will contain all needed information to extract the 
specific fuel compositions in real-time. This, how- 
ever, requires that the detector itself does not con- 
tribute any meaningful distortion to the energy 
spectrum, whether it be a non-linear response or 
a non-uniformity based on position dependence. 



II. DETECTOR DESIGNS 

The detector designs that were studied were in- 
tended to explore various methods of reducing the 
overall size of a detector while maintaining suffi- 
cient statistical acceptance to be relevant for the 
desired monitoring. The baseline fiducial volume 
of 2 m 3 (~2 tons) was chosen for all designs to pro- 
vide a sufficient event rate for any detector within 
~60 m of a reactor core. In addition, an attempt 
was made to keep the overall size of the detector 
small — all dimensions less than 3-4 m. 

The four designs, described in detail below, 
make use of the latest technological developments 
which have been implemented within the reactor 
neutrino community. All of the designs include a 
fiducial "target" region which is filled with a Gd 
doped liquid scintillator. Some of the designs in- 
clude a "gamma-catcher" which is an un-doped 
liquid scintillator region surrounding the target 
whose purpose is to completely absorb the pho- 
tons released by the n-Gd capture. Acrylic is 
used wherever a volume boundary requires optical 
transparency. Also, in all designs, 8" PMTs have 
been used and are installed within aim mineral 
oil buffer whose purpose is to reduce the singles 
background by attenuating the radioactive pho- 
tons emitted by the 40 K in the PMT glass. 



Design 1: Miniature Physics Design 

This design is a scaled down version, based on 
the state-of-the-art detectors that have been devel- 
oped for the next generation reactor based neu- 
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trino oscillation measurement at Double Chooz. 
The detector uses nested acrylic cylinders to con- 
tain the target and gamma-catcher volumes. The 
2 m 3 target has a height of 1.5 m and diame- 
ter 1.3 m while the gamma-catcher has a uniform 
thickness of 35 cm. The 1 m thick mineral oil 
buffer completely surrounds the active scintillator 
regions and is contained within a stainless steel 
vessel. The stainless steel is considered to be pol- 
ished with a reflectivity of —40%. A total of 282 
PMTs are installed on the inner wall of the stain- 
less steel vessel for an active coverage of ~15%. 

The total dimension of this detector design is 
4.2 m in height and 4 m in diameter. While this is 
a little larger than the stated goal, it is expected 
that this design will provide the optimal detector 
performance and is therefore a useful benchmark 
against which the other detector designs can be 
compared. 



Design 2: Two Sided Cylinder Design 
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This design attempts to reduce the overall foot- 
print and cost of the previous design by using re- 
flective walls to eliminate most of the PMT cover- 
age and the associated mineral oil volume. Iden- 
tical to Design 1, this design has a 2 m 3 target 
cylinder of 1.5 m height and 1.3 m diameter and a 
uniform 35 cm gamma-catcher. However, in this 
design, PMTs are only located at the top and bot- 
tom of the detector. Therefore, the outer vessel 
matches the diameter of the gamma-catcher, but 
extends an extra meter at the top and bottom to 
house the mineral oil buffer and the PMTs. The 
barrel region of the outer vessel is defined to have a 
diffuse reflective surface — approximating Tyvek — 
with reflectivity of ^90%, while the top and bot- 
tom arc defined to be non-reflective. A total of 30 
PMTs (15 each at the top and bottom) provide 
an effective covering of ^15%. The overall dimen- 
sion of this detector is 4.2 m in height and 2 m in 
diameter. 



Design 3: Two Sided Minimal Design 




In an attempt to reduce further the footprint of 
Design 2, this detector has eliminated the gamma- 
catcher entirely. To compensate for the expected 
efficiency loss at the edge of the target, the target 
volume has been expanded to match the diameter 
of the gamma-catcher in Designs 1 and 2. Thus 
the target has dimensions of 1.5 m height and 2 m 
in diameter for a total volume of 4.71 m 3 — more 
than double the previous values. This provides a 
good test of the relative benefit of increased effi- 
ciency vs. increased fiducial volume. Similar to 
Design 2, a 1 m mineral oil buffer and PMTs are 
placed at the top and bottom of the detector. As 
in Design 2, the barrel region of the outer ves- 
sel has a diffuse reflective surface consistent with 



Tyvek (90% reflectivity) and the top and bottom 
of the outer vessel are non-reflective. A total of 
30 PMTs (15 each at the top and bottom) pro- 
vide ~15% coverage. The total dimension of this 
detector is 3.5 m in height and 2.0 m in diameter. 



Design 4: Single Sided Design 




This detector is an alternate attempt to reduce 
the footprint of the detector, while maintaining 
the use of a gamma-catcher. In this design, PMTs 
are only used at the top of the detector. The tar- 
get region has the same 2 m 3 volume as Designs 
1 and 2 — 1.5 m height and 1.3 m diameter. In or- 
der to improve light collection at the single end, 
the gamma-catcher is now a conic section with a 
base diameter of 2 m, a height of 2.2 m, and a top 
diameter of 2.50 m. This provides a thickness of 
35 cm or greater around the target region. With a 
diffuse reflective surface similar to Tyvek (reflec- 
tivity ^90%) the chosen 6.5° angle of the outer 
vessel ensures that light from anywhere in the tar- 
get volume will have a path length of less than 
14 m to the plane of the PMTs. While the bot- 
tom of the detector is now also reflective, the top 
is still defined to be non-reflective. A total of 24 
PMTs are installed at the top of the 1 m mineral 
oil buffer providing ~15% active converage. This 
design has a total dimension of 3.2 m in height 
and 2.5 m in diameter. 



III. NEUTRON IDENTIFICATION 
EFFICIENCY 

The identification of an inverse beta-decay event 
is largely dependent on the ability to identify the 
signature of the neutron capture by gadolinium. 
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The n-Gd capture process releases between 3-10 
photons with a total energy of ~8 MeV. Since the 
positron energy deposition is usually between 1- 
6 MeV, it is common to define a neutron capture 
signature as any energy deposition of 6-10 MeV. 
In addition to inefficiencies which are caused by 
losses of some of the photons from the n-Gd cap- 
ture, there will also be some neutrons which will 
instead be captured on free protons in the scin- 
tillator. This n-proton capture releases a single 
photon of energy 2.2 MeV and would not satisfy 
the selection criteria. 

For this study, we generated a uniform dis- 
tribution of neutrons within the target volume. 
The neutrons were generated with a kinetic en- 
ergy of 2.5 MeV in a random direction. This ki- 
netic energy is higher than that expected from 
inverse beta-decay events, but subsequent cross- 
checks demonstrated that the results presented 
here are consistent with those of lower energy neu- 
trons. The simulation libraries described above 
were used to montecarlo the neutron thermaliza- 
tion and eventual capture. The ensuing optical 
photons from the scintillator were tracked to the 
photocathodes of the PMTs. Within the simula- 
tion, the PMTs used a radially dependent quan- 
tum efficiency which is consistent with experimen- 
tal testing. This yielded a photon hit-count per 
PMT that is roughly equivalent to a number of 
photoelectrons. This photon hit-count has been 
used for all comparisons of detector response — 
avoiding the complications of varying electronics 
and data acquisition systems. 

By using the truth information from the monte- 
carlo simulation, events were selected in which the 
neutron was stopped or captured within the tar- 
get volume. All PMT hits recorded within 100ns 
were summed together to create an effective to- 
tal charge deposition. Fig. 2] shows the photon hit 
spectrum for these events from Design 1. One can 
quite clearly see the n-Gd capture peak at around 
3200 photon hits (representing ^8 MeV) and the 
n-p capture peak at around 900 photon hits (rep- 
resenting 2.2 MeV). 

To avoid uncertainties in the quality of cali- 
brations or the variations in the linearity of re- 
sponse for the various detector designs, it was de- 
cided that the minimum threshold for n-Gd cap- 
ture identification would be defined as the number 
of photon hits located at the point | between the 
fitted peak values of the n-Gd and n-p captures. 
This value would be roughly equivalent to ~6 MeV 
and provided a very robust and uniform definition 
to be applied to all four designs. In Fig. |2 the 
photon hit spectrum for each of the four detec- 
tor designs is shown with the applicable threshold 
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FIG. 1: Total detected energy deposition for neu- 
trons captured within the target volume of detector 
Design 1. The photon hits represent the simulated 
number of photoelectrons at the PMTs. The peak at 
around 3200 photon hits represents the 8 MeV depo- 
sition due to n-Gd capture while the peak at around 
900 photon hits represents the 2.2 MeV deposition due 
to n-p capture. 



value superimposed. 

With the application of these thresholds, the 
overall neutron detection efficiency for each design 
is shown in Table [I] It was an initial surprise that 
Design 1 had such a low efficiency. After further 
investigation, it was understood that the small di- 
mension of the target and gamma-catcher allow 
many of the photons from the n-Gd capture to 
penetrate to the inactive mineral oil buffer before 
depositing their energy. In the other three designs, 
these photons are being absorbed in the outer wall 
of the vessel and some of that energy is re-emitted 
back into the active scintillator volumes. This 
quite clearly demonstrates that while Design 1 is 
optimized for larger detectors with fiducial volume 
of greater than 10 m 3 , it is not necessarily optimal 
for a more compact design. 



TABLE I: Total efficiency for selection of neutron cap- 
tures for each detector design. The sample is based on 
events in which a neutron stopped within the target 
volume. The neutron capture is identified by a num- 
ber of photon hits greater than the threshold shown in 
Fig.H 



Detector 


Efficiency 


Design 1 


51.6% 


Design 2 


80.4% 


Design 3 


50.2% 


Design 4 


83.4% 
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FIG. 2: Total detected energy deposition for neutrons captured within the target volume for each detector design, 
as labeled. The energy is shown in units of photon hits which are a montecarlo representation of photoelectrons 
from the PMTs. Also shown, for each detector design, is the threshold for selecting an event in which the neutron 



was captured on gadolinium. 
n-Gd and n-p captures. 



The threshold was chosen as the point ^ between the fitted peak values of the 



Uniformity of Neutron Identification 

To better understand the specific characteristics 
of each design, we investigated the dependence of 
the neutron capture identification efficiency on po- 
sition within the target volume. The truth infor- 
mation from the montecarlo simulation was used 
to identify the vertex position of the neutron cap- 
ture. The events were binned in the vertical (z- 
position) and radial distances from the center of 
the target volume. Using the same threshold de- 
scribed above for the definition of an accepted neu- 



tron capture event, the efficiencies, as a function 
of both z-position and radius, are plotted for all 
detector designs in Fig. |2l 

As expected, one can quite clearly see that, de- 
spite the low overall efficiency, Design 1 shows the 
most uniform response. In a similar fashion, De- 
sign 3, with the lack of any gamma-catcher, shows 
significant degradation in the efficiency as it nears 
the boundaries of the target volume. Perhaps 
more interesting is to compare the radial perfor- 
mance of Design 2 and Design 3. Their efficiencies 
are relatively similar out to the nominal radius of 
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FIG. 3: Neutron identification efficiency as a function of the position within the target volume for each detector 
design, as labeled. The truth information from the montecarlo was used to identify the Z positions (left) and the 
radial positions (right) for each neutron capture event. 
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just over 600 mm. However, in the additional tar- 
get volume of Design 3 (radius > 650 mm) the ef- 
ficiency drops off quit radically. Perhaps the most 
surprising observation is the relatively good per- 
formance of Design 4. Its z-position dependence is 
not as good as Designs 1 or 2, but for most of the 
range, its efficiency for accepting neutron capture 
events is significantly higher — 10% higher than 
Design 2 and 80% higher than Design 1. It is 
only in the bottom 25-30 cm that the efficiency 
of Design 4 drops off to values equal to those of 
Design 2. Radially, Design 4 appears to maintain 
better uniformity than all other designs with the 
exception of Design 1, only showing a declining 
efficiency for the last 10 cm. 

Of course, most of the volume of the target 
cylinder is at the larger radii where the efficiencies 
are dropping off. In fact, the efficiencies shown in 
Table [I] are under- weighted at large radii due to 
the loss of 2-3% of the generated neutrons which 
have leaked out of the fiducial volume. In a real 
detector, these events would be compensated by 
an almost equal number of events, occurring out- 
side of the fiducial volume, for which neutrons 
would be captured within the target region. To 
better evaluate the total effective volume of each 
of these designs, we assumed a uniform distribu- 
tion of neutron captures and integrated the dif- 
ferential efficiencies over the total volume of each 
design's target region with the following results: 
Design 1 = 0.93 m 3 ; Design 2 = 1.42 m 3 ; De- 
sign 3 = 2.15 m 3 ; Design 4 = 1.51 m 3 . Recall that 
the nominal target volume was 2 m 3 and that De- 
sign 3 actually had a volume of 4.72 m 3 . 



tual positron annihilation. The ensuing optical 
photons from the scintillator were tracked to the 
photocathodes of the PMTs. Within the simula- 
tion, the PMTs used a radially dependent quan- 
tum efficiency which is consistent with experimen- 
tal testing. This yielded a photon hit-count per 
PMT that is roughly equivalent to a number of 
photoelcctrons. This photon hit-count has been 
used for all comparisons of detector response — 
avoiding the complications of varying electron- 
ics and data acquisition systems. All PMT hits 
recorded within 100ns were summed together to 
create an effective total charge deposition. An ex- 
ample of the total energy deposition from these 
two samples in Design 1 is shown in Fig. 
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FIG. 4: Total detected energy deposition for positrons 
generated within the target volume of detector De- 
sign 1. The photon hits represent the simulated num- 
ber of photoelectrons at the PMTs. The positrons 
were generated with kinetic energies of 1.5 and 
3.5 MeV. 



IV. POSITRON ENERGY RESPONSE 

The positrons from a reactor induced inverse 
beta-decay event will deposit between 1-6 MeV of 
energy into the scintillator with a maximum like- 
lihood value of ^2.5 MeV. We performed some 
simple studies of the linearity of the energy re- 
sponse to positrons with kinetic energies between 
0-5 MeV and did not see any significant energy de- 
pendence. We therefore concentrated on two spe- 
cific positron samples — one with kinetic energy of 
1.5 MeV corresponding to the peak energy deposi- 
tion and one with kinetic energy of 3.5 MeV which 
would be in the high energy tail — which we used 
to study in detail the position dependence of the 
different detector designs. 

Each positron sample was generated uniformly 
within the target volumes of each detector design. 
The simulation libraries described above were used 
to montecarlo the ionization losses and the even- 



The truth information from the montecarlo sim- 
ulation was used to identify the vertex position of 
the positron event. The events were binned in the 
vertical (z-position) or radial distances from the 
center of the target volume. The recorded photon 
hit-counts were fitted to a Gaussian for each bin 
in either radial distance or z-position. These data, 
for all detector designs, are shown in Fig. [3] as a 
percent difference from the mean response of the 
detector. The errors shown represent the Gaussian 
sigma of the fit. 

The first thing to notice when looking at these 
distributions is that the responses to the 1.5 MeV 
and 3.5 MeV positrons is quite nearly identical. 
Also, the errors are consistent with the expected 
energy resolution from the liquid scintillator. As 
with the neutron capture efficiencies, Design 1 
shows the most uniform response of all detec- 
tor designs. It was a bit surprising to see the 
large variations in the radial responses for De- 
signs 2-4. After some further investigation, it 
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FIG. 5: Relative positron energy response as a function of position within the target volume for each detector 
design, as labeled. The truth information from the montecarlo was used to identify the Z positions (left) and the 
radial positions (right) for each positron event. The responses of positrons with kinetic energy of 1.5 and 3.5 MeV 
are shown simultaneously as a percent difference from the mean response of the detector. Errors represent the 
fitted Gaussian sigma of the data in each bin. 
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was understood that this structure is consistent 
with the ring-like implementation of the PMTs at 
the top and bottom of the detectors. The addi- 
tional PMTs on the barrel of Design 1 are clearly 
serving to reduce these variations. Also, as was 
expected, the detector designs which incorporate 
two-sided readouts (Designs 2-3) show more uni- 
form response relative to z-position than is seen 
in Design 4. However, given that the radial re- 
sponses of these detector designs show such large 
variations, it is not clear that the z-position unifor- 
mity provides a significant gain. In fact, it appears 
that Design 4, with a slightly larger surface area of 
PMTs at the single end, provides a more uniform 
response with respect to radial position than De- 
sign 2 — ±12% compared to ±15%. In addition, 
this appears well matched to the z-dependent re- 
sponse of Design 4 which is relatively similar. 



Simple Energy Calibration 

An attempt was made to see if the variations in 
detector response could be reduced by implement- 
ing a simple energy calibration as a function of 
vertex position. Clearly Design 1, while providing 
sufficient PMTs to reconstruct a vertex position 
with good accuracy, provides a good uniformity of 
response and was therefore not attempted. On the 
other hand, Design 4 demonstrated that readout 
from only one side was insufficient to consistently 
reconstruct a vertex position that would be useful 
for correcting the energy response. 

For Designs 2 and 3, we attempted various 
methods to reconstruct the radial position. How- 
ever, the reflectivity of the side walls made this 
unreasonably difficult. We did find, however, that 
a simple charge balance between the top and bot- 
tom PMTs provided a very good correlation to the 
true z-position. We define the measurable quan- 
tity z-balance (Zb a i) in the following way: 



Zbal 



Top — Bottom 
Top + Bottom ' 



(1) 



where Top and Bottom refer to the total num- 
ber of photon hits recorded by the PMTs at the 
top and bottom of the detector, respectively. The 
correlation of this value to the true z-position for 
Designs 2 and 3 can be seen in Fig. Eland the rel- 
ative positron energy response (similar to Fig. [SJ 
is plotted relative to this variable in Fig. [7| Us- 
ing the data from Fig. [7] a multiplicative calibra- 
tion constant for each bin in Z^ai was constructed 
to give the mean response. Rather than creating 
an energy dependent calibration, the relative re- 
sponses of the 1.5 and 3.5 MeV positron samples 
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FIG. 6: Correlation between the measurable z- 
balance, as defined in Eq. and the true z-position 
from the montecarlo for detector Designs 2 and 3. 
While the correlation is relatively linear for both de- 
signs, one can see the non-linearity of Design 3 at the 
extremes due to the energy leakage in the absence of 
a gamma-catcher 



were averaged to provide a single calibration con- 
stant for each bin. This procedure was repeated 
for both detector designs and the calibration con- 
stants were applied event-by-event. 

As expected, after calibration the uniformity of 
the response for both designs is improved (See 
Fig. |HJ). Especially in the region of Zt a i between 
-.3 and .3, the relative differences are kept within 
±1% — consistent with the uncalibrated perfor- 
mance demonstrated in Design 1. One can also 
note the larger error for Design 3 which arises from 
the large response variation with radial position. 



V. CONCLUSION 

This work was an attempt to better understand 
the design features that would affect the response 
of a compact liquid scintillator detector to anti- 
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FIG. 7: Relative positron energy response as a func- 
tion of the z-balance as defined in Eq. for detector 
Designs 2 and 3, as labeled. While the gross structure 
of the position dependence is similar to that shown in 
Fig. |S] the magnitude of the variation is slightly less. 

neutrinos from nuclear power reactors. There are 
some interesting trade-offs that can be identified. 
For example, Designs 2 and 4 offer an interesting 
compromise in efficiency and uniformity while De- 
sign 1 shows that very good uniformity of response 
can be achieved at the cost of a larger total detec- 
tor size and an efficiency for event selection of only 
50%. On the other hand, Design 3 demonstrates 
the maximum effective fiducial volume for the min- 
imum total detector size at the cost of energy dis- 
tortions as large as 40%. The optimal design will 
depend on the goals and restrictions of the specific 
installation. 

Clearly improvement can be made to any and 
all of these designs. It was perhaps interesting to 
note that the expected uniformity from the two- 
sided designs was hampered by the large radial 
variations in energy response. Some more clever 
design of the PMT distribution might be able to 



improve this. Similarly, one might consider some 
judiciously placed additional PMTs in Design 4 to 
allow for some minimal reconstruction and cali- 
bration. 

In general, however, it appears that a compact 
neutrino detector with overall dimensions on the 
order of 3 m can be built with reasonably good 
performance. Given that the energy resolution of 
the scintillator is usually 7-10%, the 15% variation 
in energy response is not extreme. It is the hope 
of the authors that this work will prove useful to 
those working on nuclear reactor monitoring and 
can provide fruitful directions for future detector 
development. 
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FIG. 8: Relative positron energy response as a func- 
tion of the z-balance as defined in Eq. for detector 
Designs 2 and 3 after calibrations were applied. Multi- 
plicative calibration factors were derived for each bin 
by averaging the responses of the 1.5 and 3.5 MeV 
positrons shown Fig.Q These calibration factors were 
then applied event-by-event. A factor of two reduc- 
tion in the non-uniformity of response was achieved 
for both detector designs. 
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